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Abstract 

A status report on the theory of radiative rare B decays in the standard 
model is presented, with focus on inclusive decays B — > X^ s ^ 7 and exclusive 
decays B — > (K*,u,p)j. CP asymmetries are also briefly discussed. 

1 Introduction 

The radiative rare B decays have received a lot of attention in the recent years. Since 
they are forbidden at tree level in the standard model, they can occur only as loop 
effects and therefore give information on masses and couplings of virtual particles 
running in the loops, like the W or the t quark. Radiative rare B decays test QCD 
corrections and provide a searching ground for non standard physics and CP violating 
asymmetries. 

Here a status report on the theory of radiative rare B decays in the standard model 
is presented. The effects of the strong interactions on the weak radiative B decays 
are studied in the framework of the effective hamiltonian. We discuss separately the 
inclusive decays B — > X( Sjd ) 7 and the exclusive decays B — > (K*, u>, p) 7, and then 
briefly comment on CP violation asymmetries. 

1 Invited talk at European Physical Society Conference on High Energy Physics, Brussels, July 
1995; to be published in the proceedings. 
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2 Inclusive Decays B —> Xu^pf 
2.1 Rate 

The inclusive decay B — > X s 7 is described at the partonic level by the weak decay 
b — ► S7, corrected for short- distance QCD effects. Support for the spectator model in 
inclusive decays comes from the l/ra& expansion JT|. The perturbative QCD correc- 
tions are important in this decay, enhancing the rate by 2-3 times, which makes the 
theoretical prediction compatible with the experimental rates within the errors. 

The perturbative QCD corrections introduce large logarithms of the form a™ (/i) log m (/i/M) 
(m < n), where a s is the strong coupling, M is a large scale (M = m t or My/) and 
H is the renormalization scale. By using renormalization group equations, the large 
logarithms are resummed order by order and the coefficients of the effective hamil- 
tonian can be calculated at the relevant scale for B decays fi ~ 0(m.&) @. Although 
the first analyses date back to many years ago, the first fully correct calculation at 
leading order (LO) of the anomalous dimension matrix has been obtained only in 
19930 and confirmed last yearQ. The main problem has been the evaluation of the 
two loop diagrams, that mix the operators (Qi-.-Qe) with the operators (Q7, Q 8 ) (for 
the definition of Qi-..Qs see e.g. 0). The effect of these diagrams has been found too 
large to be ignored. It should not be surprising that two loop diagrams are already 
present at the LO in QCD corrections, since this weak decay starts at one loop at 
order a°. 

The next-to-leading order (NLO) calculation has been only partially completed. 
The two-loop mixing in the sector (Qi...Qe) has been calculated 0, as well as the 
two-loop mixing in the sector (Q7,<5s)0- Gluon corrections to the matrix elements 
of magnetic penguin operators have also been calculated 0, §, [|. The 0(a s ) cor- 
rections to Ci(Mw) and C$(Mw) have been considered in ref.HlQj. The three loop 
diagrams that mix the operators (Qi-.-Qe) with the operators (Q7, Q$) are still to be 
calculated; as seen at LO, their contribution may be relevant and estimates based on 



the incomplete NLO calculation must be handled with care. At LO||ll||, it has been 
estimated 

Br(B -> X sl ) TH = (2.8 ± 0.8) x 1(T 4 , (1) 

assuming |V^|/|Vd,| = 1, as suggested by the unitarity of the Cabibbo Kobayashi 
Maskawa matrix (CKM). After the inclusion of 0(a s ) virtual and bremsstrahlung 
corrections and taking into account the scale dependence of the running quark masses, 



the prediction is [12] 



Br(B -> X sl ) TH = (2.55 ± 1.28) x 1(T 4 . (2) 
The error is dominated by the uncertainty in the choice of the renormalization 



scaleP, 11, 12|. An attempt to apply scale fixing methods to this decay has been 
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done [13 1 . However, the last word will belong to the complete NLO calculation, that 
should considerably reduce theoretical uncertainties at LO. Within the large errors, 
the prediction is in agreement with the inclusive CLEO data flTJ 



Br(B -> X sl ) EXP = (2.32 ± 0.57 ± 0.35) x 10" 4 . (3) 
By factorizing the CKM parameters, the amplitude can be written as 

A = v u A u + v c A c + v t A t , (4) 

where v u = V* s V u b, v c = V* s V c b and v t = V t * s V t b- Since v u is negligible with respect to v c 
and v t (\v u \ ~ 0(A 5 ) and \v c \ ~ \v t \ ~ 0(A 2 ), in the parameterization of Wolfenstein) , 
it follows that v c ~ — Vt by the unitarity of the CKM matrix. The amplitude is thus 
proportional to v t and one can use this decay normalized to the semileptonic decay 



to estimate \Vt s \/\ Vcb\ W%\ 



\Vu\l\Vcb\ = 1.10 ±0.43. (5) 

Similarly, one could use b — > d'-f to extract information on \Vtd\- This decay has 
not been detected. The expected branching ratio is approximately O(10~ 5 )[|l5|. Even 
though statistically the inclusive decay is in the reach of future CLEO /B factories, it is 
difficult to be observed, due to the large background of b — > 57. In b — > dj, the CKM 
factors in the effective hamiltonian have comparable size |i> u | ~ \v c \ ~ \v t \ ~ 0(A 3 ); 
the proportionality to \V t d\ could thus be jeopardized by contributions coming from 
the c and u loops. At LO, however, there are no large contributions of the type 
a s log(m^/rn|)fT§. 



Recently, the possibility of non perturbative long distance effects on the rate 
through resonant intermediate states has been suggested. These effects have been 
estimated by using the vector meson dominance hypothesis (VMD) or non relativistic 
quark models[|l6|, [17], ^8[. In the estimates by VMD, the radiative transition b — ► 
(s, rf)7 is modelled by a sum over the processes b — > (s,d)V*, where V* is a virtual 
vector meson (tp and its excited states, p, and a>), followed by the conversion V* — > 7. 
In this picture, it is assumed that only the transverse part of b —>■ 57 couples to the 
photon, in order to preserve gauge invariance. A crucial problem is that the effective 
couplings, that are measured at the mass of the resonance in the intermediate process, 
must be scaled to q 2 = 0, since the final photon is on shell. This scaling may lead to 
a strong suppression that reflects in small long distance contributions, less than 10% 



with respect to the short distance rate [16, 17]. The size of this suppression awaits 



further work; if absent, the rate becomes much bigger fTq] . Similar results follow in 



the b — ► g?7 case|16, 17 



2.2 Photon energy spectrum 

In the two body decay b — >• 57 the photon energy is fixed. A non-trivial photon energy 
spectrum is obtained by 
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including perturbative emission of hard gluons, such as b — > s'jg [15 



taking into account the non perturbative motion of the b quark inside the 
meson 



12, 19, 20 



The amplitude for the decay b — > s^g suffers from singularities in the limit of 
soft gluons or photons (E 7 — > E™ ax or E^ — > 0, respectively). These singularities are 
cancelled in the photon energy spectrum if one also takes into account the virtual 
corrections to b —>■ sj and b — > gj, order by order in perturbation theory Near 
the endpoint regions, the spectrum can be improved by resumming at all orders the 
leading (infrared) logarithms fi~2|, p0| . The region E 1 — > E™ ax deserves particular 



attention, since it is the region that contributes mostly to the rate. Note that in the 
limit m s — > also collinear singularities come into play|[21||. 



In order to implement the 5-meson bound state effects on the photon energy 
spectrum, one can use a specific wave function model p2| , where the 5-meson consists 



of a bound state of a b quark and a spectator quark of mass m q . The b quark is given 
a momentum having a Gaussian distribution, centered around zero, whose width is 
determined by a parameter pp. Both parameters of the model, pp and m q , can be 
fitted by the CLEO photon energy spectrum fT2| . 



Another approach to the B- meson bound state effects is based on QCD and l/m& 
expansion |19"|, pO|j . The spectrum is expressed in terms of a universal distribution 
function whose moments are related to local quark operators of increasing dimensions. 
This function depends on the final quark mass (e.g., it is different for B — > X s ~f and 
B — > X c eu). A few free parameters have to be determined by matching with the 
experimental data. 

In inclusive b (s, d)~f decays, the two approaches are compatible through the 
leading order in 1/mJ 



3 Exclusive Decays B — > Vj 



The matrix element of the effective hamiltonian gives the so called short distance 
(SD) contribution to the amplitude for the exclusive decays. For the B — > Vj decay 
(where V = K*,p,u), the SD amplitude A is proportional to the matrix element of 
7 



< V(r), k)\sa 



1+75 



q v b\B{p) >= 2T^ v (q 2 ) e^ Va (k)p p k a + 



(q 2 ) mk){m% - m 2 v ) - ( V (k) ■ q){p + kf] + 



iT, 



B^V 



(Q 2 ) (V(k) ■ q) 



m 



B 



{p + ky 



(6) 
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where T^ v (q 2 ) , T^^ v (q 2 ) and T^^ v (q 2 ) are real form factors, rj is the V polariza- 
tion vector and q = p — k is the photon momentum. One can show that when the 
photon is on-shell T-f~ >v (0) = T^^ v (0) and T^ v does not contribute to A; the rate 
depends thus on one form factor only. 

Among exclusive decays, the B — > K*^ has received an increasing attention and 
the relative form factor has been calculated in a large number of papers (see e.g. 
and references therein), employing several models: HQEFT, quark models, QCD sum 
rules, lattice. When comparing experiment and prediction, it is convenient to use the 
ratio Rk* = T(B — > K*j)/T(b — > 57), that is largely independent from many the- 
oretical uncertainties, like renormalization scale, unknown perturbative higher order 
corrections, etc. Experimentally, Rk* = (19 ± 9)%[p5|, |14|| . Recently, QCD sum 
rule models have given results in the ball park of the experimental data[^, p7| , e.g. 
Rk* = (16 ± 5)%P7[ by using the light cone sum rules. 



In lattice QCD, the point of physical interest is not directly accessible at present 
and one has to make an ansatz about extrapolating the results to q 2 — > and to the 
physical b quark mass. In particular, the functional dependence of on Q 2 is critical 
for the results, since a pole-like behaviour for T^^ K * (q 2 ) results in Rk* ~ (4 — 13)%, 
while a constant behaviour for T^~* K (q 2 ) results in an appreciably higher rate. There 
are four groups working on this matrix element: BHS[^|, UKQCD|29[, LANLpU| and 



APE|3l|; we limit here to report their results for Rk* in the cases of T^ K {q 2 ) pole 



dominance (first number) and T^^ K * (q 2 ) constant (second number): 

BHS: (6.0 ±1.2 ±3.4)%, (7) 

UKQCD : (13 + 14 - 10)% (35 ±4 -2)%, (8) 

LANL: (4-5)% (27 ±3)%, (9) 

APE: (5 ±2)% (31 ±12)%. (10) 



The exclusive decay B — >• p7 has not been seen yet, but it is likely to be seen at 
future CLEO/B factories. It can be compared to B — > K*^ to extract information 
on ivy/ivy. Although the form factors of B — > K*^ and B — > p7 decays are model 
dependent, their ratio should be more reliable, being determined by 577(3) symmetry 
considerations. Therefore, in the limit of SD contributions only and assuming that 
the top quark loop dominates the ratio, | V^|/| Vt s | can be estimated by 



T(B M ^ P1 ) _\V t 



td\ 



12 



r(5 (M) - #* 7 ) \v t 



ts 



12 



eo, (11) 



where ^ is the squared ratio of the form factors and Q is a phase space factor. The 
estimate is[p2] 

< (0.64- 0.76), (12) 



t.s 
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where the range reflects the model dependence. 



Until now we have assumed that only the matrix element of the magnetic moment 
operator Oj contributes to the rate, neglecting the possibility of other long distance 
effects. In particular, beyond the spectator model approximation, the diagrams where 
the b quark annihilates the spectator quark by weak interaction have also to be 
considered. Because of the CKM matrix elements, this annihilation mechanism is 
negligible for B — > K*^, but it may be important for other exclusive radiative decays 
like B — > /7y. The presence of these effects may invalidate the relations ([fl]) and fll2[) . 

Due to color suppression, the contributions of weak annihilation diagrams in neu- 
tral B decays B° — > uj) 7 are generally believed to be smaller than in the B^ 
case. Non spectator diagrams for B ± — > p ± 7 have been evaluated in a constituent 
quark model and found to change the decay rate by a factor of 0.7 — 2.5|$3| with 
respect to the SD contribution. Estimates based on QCD sum rules find that the 
contribution of the weak annihilation diagrams modifies the SD rate of B^ — > p ± + 7 
up to ±20% I 



The analysis of long distance contributions by VMD in exclusive decays B — > 
(K*,p)^ presents many theoretical uncertainties, like e.g. the role played by the 
spectator quark. These uncertainties reflect in a wide range of results for the VMD 
amplitude; it has been estimated to be from 5% to 50% of the SD amplitude [|18|, |3~6| , 



We stress that in exclusive decays particular care must be exercised to avoid 
possible double counting among long distance effects. 



4 CP Violation 

CP violation in B radiative penguins decays may occur as interference among loop 
diagrams involving the u,cort virtual quarks. Gluon exchange provides the necessary 
strong phase shifts between these diagrams. By using the unitarity of CKM matrix 
(yTs^tb = —Vu S Vub — Ks^cb)) we can write the amplitude in the form 

A — v u A u + v c A c , (13) 
where v u = V* s V u b and v c = V* s V c b. The CP violating asymmetry is 

rcg -> J) - t(b -. /) 
acF ~ rW^TJTTW^Ty ( ' 

where B — > f and B — > f are CP conjugate processes. The asymmetry fll4]) can be 
written as 

= -4Im(^ c *) lm(A u A* c ) 
a ° P \v u A u + v c A c \ 2 + \v*A u + v*A c \ 2 ' 
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The amplitude Ai has a strong phase that does not change sign in the transformation 
Ai — > Ai, in contrast to the weak phase due to the Cabibbo Kobayashi Maskawa 
factors. At SD, the amplitude A is the matrix element of the effective hamiltonian 
between the initial and final state. At LO, only O7 contributes to the decay, the 
amplitude is real and the asymmetry is zero; one needs to go at least at order 0(a s ) 
in the matrix elements to create the phase difference between A u and A c . This 
asymmetry has been estimated |38] to be of the order (0.1 — 1)% for b — * 57 and 



'1 — 10)% for b — > dj. For exclusive modes, a C p is typically 1% for B — > K*j and 
15% for B — > pry^RQ. It is evident that the observation of a large asymmetry in 
b — > S7 mediated processes would provide by itself strong evidence of physics beyond 
the SM. 

The asymmetry (0) has been estimated in a constituent quark model, including 
the contributions due to non spectator diagrams, for B ± — > p ± 7; it has been found 
that it can be sizable, possibly as large as 30%|3"3"|. 



The prospects for CP violation in B decay modes that are dominated by at least 
two interfering resonances have also been investigated [fiOfl . It has been estimated that, 
in decays like B — > Kn^, B — > K*n r y ) B — > Kfry and -B 1 * 1 — > 7T ± 7r ± 7r ± 7, CP violating 
distributions may be observed in a sample of about 10 8 — 10 9 B ± mesons [|0 . 
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